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Abstract. Time- varying gravitomagnetic fields are considered within the linear post- 
Newtonian approach to general relativity. A simple model is developed in which the 
gravitomagnetic field of a localized mass-energy current varies linearly with time. The 
implications of this temporal variation of the source for the precession of test gyroscopes 
and the motion of null rays are briefly discussed. 
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1. Introduction 

Astronomical sources, such as neutron stars, generally rotate, but their rotation rates 
are seldom uniform. For instance, the Earth's rate of rotation decreases very slowly 
due mainly to tidal friction; in fact, in this process the Earth loses angular momentum, 
which is then transferred to the orbital motion of the Moon around the Earth. The 
temporal variation of the proper angular momentum of a massive body in turn generates 
a time- varying gravitomagnetic field in accordance with the general theory of relativity. 
The purpose of this paper is to study gravitomagnetic fields generated by the intrinsic 
temporal variability of sources; this issue was ignored in a previous work on general 
gravitomagnetic fields of spinning masses pQ. Only a beginning is made here, since the 
detailed study of the gravitomagnetic field resulting from a realistic model of a variable 
source is beyond the scope of this first analysis. 

The linear approximation scheme of general relativity is adopted throughout this 
paper. Thus the spacetime metric tensor g^ v is expressed as g^ = T]^ + h^, where 
r]^ is the Minkowski metric tensor with signature +2 and h^ u is a first-order ( "weak" ) 
perturbation. The choice of the background global inertial coordinates x^ = (ct, r) 
is not unique; under a slight transformation of coordinates x M t— > x^ 1 — e M , we have 

t— > h^ v + e^ u + e„ iM , so that the gravitational potentials are gauge- dependent 
quantities. However, the Riemann, Ricci and Einstein curvature tensors are gauge 
invariant. In particular the Einstein tensor may be written as 

G,*, = - l -Uh^ (1) 

where the trace-reversed potentials h^ v are given by h^ u = h^ u — \rj^ u h. Here h = tr(/i Ml/ ), 
□ = rf^d a dfi and the transverse gauge condition 

W\ v = (2) 

has been imposed. The gravitational field equations, G^ u = ( 87rG/ c 4 )T^ ul can then be 
expressed as 

□V = -^-T^. (3) 

The solution of this linear inhomogeneous equation is a superposition of the retarded 
solution 

c 4 J |r — r | 

plus the general solution of the homogeneous wave equation. The latter is ignored in 
this analysis and we simply focus on Eq. (jlj) in what follows. Given a consistent model 
for the source T^ u , d v T^ u = 0, h^ u is thus determined and then 

hfj, v = h^ v — —f]f_iuh, (5) 



where h = tr(h 



flU } 
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It is now necessary to introduce the slow-motion assumption that is needed in our 
linear post-Newtonian approach to gravitoelectromagnetism (GEM). We require that 
all motions within the source take place with speeds that are very small compared to c. 
Furthermore, all contributions to the metric tensor of 0(c~ 4 ) are neglected. It follows 
from Eq. (Tj0) that for the sources under consideration here, /too = 4$/c 2 , hoi = —2Ai/c 2 
and hij = 0(c~ 4 ). In fact, $(t,r) is the gravitoelectric potential and A(t, r) is the 
gravitomagnetic vector potential such that 

!f + v. A = o, (a, 

and the metric is given by 

ds 2 = -c 2 (l - ^-J dt 2 — ^(A • dr)dt + (\ + dijd^dx 3 . (7) 

Our definition of GEM potentials is based on a special convention that preserves 
the electromagnetic analogy as much as possible [2]. Moreover, it is useful to define 
= (2$, A), so that h ^ = 2A^/c 2 , and the GEM Faraday tensor is then given 
by Ffjp = A V41 — A^^. It follows from these defintions and Eq. ([3]) that the Maxwell 
equations for the GEM fields are satisfied [2]; that is, \— > (2E, B), F^p] = and 

F^ = ^f, (8) 

where is the mass-energy current such that 

cf = (T 00 ,T 0i ). (9) 

Thus in the GEM approach adopted in this paper, T*- 7 does not explicitly appear; 
however, its existence is implicitly assumed such that the four dynamical equations 

T^ v = (10) 

are always satisfied. The linear perturbation approach to GEM has been reviewed in [3]. 
The compact notation of Eqs. (|S]) and has been employed for simplicity, so that the 
tensorial character of these equations is purely formal; that is, the background global 
inertial frame is essentially fixed in this linear GEM treatment. 

To study time-dependent gravitomagnetic fields caused by the intrinsic variability 
of the source, it would be necessary to have a consistent model of a source with time- 
varying proper angular momentum. This is rather complicated in general; however, a 
drastic simplification occurs within the linear GEM framework via a simple toy model 
that is the subject of the next section. 



2. Toy model 

It is an immediate consequence of Eqs. (pO) and (jSJ) that if T^ v is independent of time, 
then h^y is only a function of spatial variables and hence the spacetime of a time- 
independent source is in general stationary. Let T* represent the energy-momentum 
tensor of a time-independent source and h* be the corresponding stationary solution 
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of Eq. (j3J); then, the time- independent GEM potentials are $*(r) and ^4*(r), while 
h*j = 0(c -4 ). A detailed treatment of such configurations is contained in [I]. For time- 
varying GEM potentials, a simple ansatz would involve a separation of temporal and 
spatial variables. Under what conditions can 

r) = <p(t)$*(r), A(t, r) = f(t)A*(r) (11) 

represent acceptable GEM potentials? Equation ([6]) implies that dip/dt = 0, while we 
find from direct substitution of Eq. pip in Eq. ([3]) that d 2 f/dt 2 = and 

T 00 = ^T* , T 0l = f(t)T* t . (12) 

It follows that our ansatz ( iTTjl is consistent with GEM equations provided that is 
constant and /(£) is a linear function of time such that 2|A| << c 2 ; moreover, the 
corresponding Tij,djT- = (df/cdt)^, is expected to lead to hij = 0(c -4 ) via Eq. (TJJ. 

The mass-energy current (Q for the case of Eq. ffl2|) is such that j, cj % = T° l , is 
localized and divergenceless; therefore, its volume integral vanishes. This implies that 
Eq. (fTTj) is consistent with the retarded solution (j3J), since f(t) is just a linear function 
of time. In addition, stresses T^- exist such that the source has proper dynamics in 
accordance with Eq. ( flOl) ; however, the GEM metric (JTj) depends only on $ and A and 
hence the nature of the stresses is of no consequence here. Specifically, these stresses 
may turn out to be physically unrealistic; however, this is mitigated by the circumstance 
that the source leads to a simple analytic solution that can be consistently employed 
within the linear GEM framework. This is illustrated by the special example that is 
presented in the following section. 

With ip constant and f(t) linear in time, the GEM fields are 

1 df 

E(r) = -y, W - --LA*, B(t, r) = f(t)V x A*, (13) 

so that E is independent of time and only B is linearly dependent upon time. Thus the 
gravitational analogue of the displacement current vanishes in this case, but our simple 
model is ideally suited to help elucidate the gravitational analogue of the Faraday law 
of induction [5]. Moreover, the linear dependence of the source upon time is consistent 
with the absence of gravitational radiation in this case. 

Many significant relativistic effects of astrophysical interest fall within the linear 
GEM scheme; therefore, our simple model can provide useful results. In practice, the 
slow temporal variation of the angular momentum of the source may not be linear; 
however, our linear model could furnish a first approximation. Some applications of 
our model are discussed in the rest of this paper using a rather special nonstationary 
spacetime. 

3. Special nonstationary spacetime 

Far from any localized stationary source of mass M and angular momentum J*, the 
spacetime metric can be expressed in the form (JTj), where the dominant contributions 
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to the GEM potentials are given by 



GM GJ* J x r 
$ = , A* = — (14) 

r c r 6 

for r >> GM/c 2 and r >> J*/(Mc). Here J is a fixed unit angular momentum vector 
of the source. With these potentials, Eq. ([7]) is equivalent to the linear form of the Kerr 
metric in isotropic Cartesian coordinates. 

It follows from the results of the previous section that within the GEM framework, 
metric (J7J) with 

J= GM A= GJWJx i _ 

J(t) = f(t)J* ' (16) 

represents the exterior field of a localized source of mass M and angular momentem J(t) 
that varies linearly with time along the fixed axis J. Substituting Eqs. ffT5]) and ffTBT) in 
Eq. ([3D, one can find the effective source for these potentials; indeed, 

Too = Mc 2 5(r), T 0i = ±c[J(t) x V] 4 5(r). (17) 

Thus the mass-energy current is singular in this case and essentially confined to the 
origin of spatial coordinates. On the other hand, Eq. ffTUl) implies that 



13 8tt 



(j x V)iVj- + (J x V)jV<- 



where J = dJ/dt is a constant vector. The stresses (ITS]) are time-independent and 
extend throughout space, but fall off to zero as r~ 3 for r — ► oo. As emphasized before, 
the influence of the stresses on the spacetime metric is of 0(c~ 4 ) and is therefore ignored 
in the linear perturbation approach to GEM. 

The localized mass-energy current that generates exterior GEM potentials (1151) 
in an astrophysical situation could be considered as a variant of Newton's experiment 
with the rotating bucket of water but on a vastly different scale. Imagine an initially 
stationary configuration with angular momentum Jj. At some initial time tj, the system 
undergoes a linear spin-up (or spin-down) that ends at tp and the system then returns 
to a stationary state for t > tp with angular momentum J p. The total mass-energy of 
the configuration remains the same throughout. Thus Eq. (fT51 holds for all time and J 
is constant except for t G [i^ip], where 

J=Ji+'^—+ L (t-ti). (19) 
tp — tj 

Moreover, J(t) is piecewise constant and hence 

J = T^—riKt - ti) - 6(t - t F )]. (20) 
tp — tj 
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This implies, via Eq. Q, that there are additional spatially extended but instantaneous 
currents at tj and tp responsible for respectively turning on and off the temporal 
variation of the source. That is, Eq. ( TlT|) holds except that T 0i is modified at tj and tp, 

+ ic(Jx V)itf(r). (21) 

It is interesting to explore the influence of time-varying gravitomagnetic fields on 
some standard consequences of general relativity. In the following subsections, we 
employ the nonstationary linearized Kerr spacetime to discuss some aspects of spin 
precession and the motion of null rays; the motion of test particles will be treated in 
[5] . For the sake of simplicity, we limit our considerations to temporal intervals that are 
within and are such that 2|A| << c 2 . 




3.1. Gyroscope precession 

The purpose of this subsection is to provide an estimate of the influence of intrinsic 
variability of the source on gyro precession. A detailed treatment of gyro motion in 
the special nonstationary gravitational field under consideration here on the basis of 
the Mathisson-Papapetrou-Dixon equations is clearly beyond the scope of this work [6]. 
Instead, we concentrate on a free test gyro held fixed outside the source with potentials 
given in Eq. ffl5l) . It follows from a simple extension of previous results j6] that to linear 
order in the spin, the gyro precesses with frequency B/c, where the gravitomagnetic 
field is given by 

B = ^P[3(J-r)r-J]. (22) 

The resulting motion can be simply described via introduction of a new temporal 
coordinate r given by dr = f(t)dt. 

For a test gyro in orbit about the source, the motion of the gyro spin is more 
complex; nevertheless, it is possible to estimate the influence of a variable mass- 
energy current on the gyro spin. For the Earth, the secular increase in the period 
of proper rotation P e is approximately linear in time and can be characterized by 
~ 3 x 10 -10 per year. This implies that the influence of the secular decrease of the 
Earth's angular momentum on the GP-B experiment [7J is entirely negligible. Indeed, 
the magnitude of the corresponding net decrease in the strength of the gravitomagnetic 
field of the Earth over the course of the GP-B experiment is about eight orders of 
magnitude smaller than the projected experimental sensitivity of the GP-B [TJ. 

The spin-curvature force on the test gyro is given by 

F a = -\R a P ,yS^\ (23) 

where u 13 is the four- velocity of the gyro [6] . It follows that in the linear GEM framework 
and to first order in spin, F° = and 

C F = -V(S ■ B). (24) 
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Thus the classical spin-gravity coupling is recovered, except that now B is linearly 
dependent upon time. We expect that the same holds for intrinsic spin with time- 
dependent Hamiltonian H = S • B/c; time-dependent spin- rotation coupling has been 
considered in [8]. It is straightforward to verify this coupling for photon spin via the 
rotation of plane of polarization of electromagnetic radiation propagating in the exterior 
field of the source. The general scheme given in [TJ implies that the result is 

9 = 1 J B • c2r, (25) 

where the integration is carried out along the path of the null ray. 

Imagine an electromagnetic ray that propagates along the axis of rotation of the 
source. During the time that the ray travels from r\ to r 2 , the angular momentum of the 
source changes linearly from J\ to J 2 and the plane of polarization of the electromagnetic 
wave rotates by an angle 9, where 

d G {r 2 -n) (Ji + JA (26) 
c i rir2 y ri r2 j 

The standard Skrotskii result is recovered for J\ = J 2 pQ. 
3.2. Time delay 

We next examine the gravitomagnetic time delay of null geodesic rays in our model. 
The results may be of astrophysical interest in pulsar timing, for instance, as pulsars 
generally lose angular momentum due to external electromagnetic braking torques. 

Within the linear GEM framework, the time delay Aq of a null ray propagating 
from Pi : (cii,ri) to P 2 : (ct 2 ,r 2 ) is 

t 2 -t 1 = -\r 2 -r 1 \+A G , (27) 



c 



where [9] 



'Pi 

Here k a = dx a /d( = (1, k) is the constant tangent vector of the unperturbed null 
ray and the integration is performed along the straight line from Pi to P 2 . Thus 
Ag = Age + Agm, where Age is the familiar Shapiro time delay and 



2 



*GM 



C 



P2 

A ■ dr (29) 

Pi 



is the gravitomagnetic time delay [9]. 

It is interesting to work out A G m for the exterior field given by Eq. ([151) . The result 

is 



2GJ- (h x r 2 ) f J x , J 2 



c 4 1 + f 1 • r 2 V r i r 2 



- + - , (30) 
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where Eqs. 0A.30 and OA .40 of the Appendix have been used. When the ray propagates 
from Pi to P 2 , the angular momentum of the source varies linearly from J\ to J 2 ; for 
J\ = J 2 however, Eq. ( 1301) reduces to formula (16) of [9]. 

The time delay of null rays in the presence of a cosmological constant has been 
similarly calculated in [10J. This treatment has been erroneously criticized in a recent 
paper [11]. Contrary to the claim in [11] , the argument in [10], as is evident from the 
context, was simply based on the mathematical fact that for a null ray in Schwarzschild- 
de Sitter spacetime in standard coordinates, the geometric path of the ray in space does 
not suffer an additional deflection away from a straight line due to the presence of a 
cosmological constant — see, for instance, Eq. (9) of |llj . 

3.3. Deflection 

Finally, for a null geodesic ray with tangent vector that travels from Pi to P2, 
K tl (P 2 ) - K^Pi) = where 

E" = - f\%k a ^dC, (31) 
Jpi 

follows from the integration of the null geodesic equation in the linear GEM scheme. 
Using the connection coefficients given in [3] for metric ([7]) with potentials ffl5l) . it is 
straightforward to compute the integral in Eq. ( !3~TT) along the unperturbed path of the 
ray. Thus S M = (E°, S) is given by 

£° = -^(QJo+Xi) 

cr 

an 1 

+— (n x k) ■ [j^QCo + A) + - + £ 2 )}, (32) 

& c 

2GM - 2G 

S = — (r x - Qk)X - -^r(ri x J 2 - r 2 x J^Jq 

an 1 

+— n x k [Ji • ( ri £ + k£i) + -J ■ (ri£j + k£ 2 )]. (33) 
Here Q = ri • k, L = |r 2 — ri| and 

k=^(r 2 -n), j = j{J 2 - Jx)J. (34) 

The integrals X , X 1; £ , £1 and £ 2 are given in the Appendix. It is interesting to note 
that 

If t 

s • k = -j-n • (k x J)J . (35) 

The result of this calculation may be of interest in connection with the significant 
advances that are expected to occur in microarcsecond astrometry in the near future [12J. 



4. Discussion 



The gravitational field equations have solutions in which the mass-energy and angular 
momentum of the source vary with time as a consequence of emission or absorption of 
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radiation. However, the present paper is devoted to nonradiative situations involving the 
slow temporal variation of the angular momentum of the source. A method is presented 
to estimate the influence of such variability on gravitomagnetic phenomena. 

The main result of this work is based on the observation that within the linear 
GEM framework it is possible to transform a localized stationary configuration to 
one with a gravitomagnetic field that varies linearly with time. Starting with 
the special nonstationary spacetime corresponding to the linearized Kerr metric in 
isotropic Cartesian coordinates, the consequences of the linear temporal variation of 
the gravitomagnetic vector potential are explored in some situations of physical interest 
such as the GP-B experiment, the rotation of plane of polarization of electromagnetic 
waves and the time delay as well as deflection of electromagnetic signals. However, a 
proper discussion of the methods needed to deduce observable effects from the results 
presented here is beyond the scope of this work. 



Appendix 

The purpose of this appendix is to evaluate the integrals 

for n = 0, 1 and m = 0, 1, 2. Here 

r(C) = r 1+ kC, (A.2) 

while k and L have been defined in section [3] These integrals can be computed using 
formulas 2.263 and 2.264 given respectively on pages 82 and 83 of [T5] . We find that 

J = ik-(f 2 -fi), X 1 = ^(l-f 1 -r 2 ), (A.3) 
o o 

where 5 = r\ — Q 2 , Q = k ■ r x as before, and hence 

5L 2 = r\r\\\ - (f x • f 2 ) 2 ]. (A.4) 

Furthermore, 

A = - Q£ - \ (\ - ^ ) , (A.6) 



A = ~r?A, - J<2A - (A.7) 



2r. 
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